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Weyl semimetal is a new topological state of matter, characterized by the presence of nondegen- 
erate band-touching nodes, separated in momentum space, in its bandstructure. Here we discuss 
a particular realization of a Weyl semimetal: a superlattice heterostructure, made of alternating 
layers of topological insulator (TI) and normal insulator (NI) material, introduced by one of us 
before. The Weyl node splitting is achieved most easily in this system by breaking time-reversal 
(TR) symmetry, for example by magnetic doping. If, however, spatial inversion (I) symmetry re- 
mains, the Weyl nodes will occur at the same energy, making it possible to align the Fermi energy 
simultaneously with both nodes. The goal of this work is to explore the consequences of breaking 
the I symmetry in this system. We demonstrate that, while this generally moves the Weyl nodes to 
different energies, thus eliminating nodal semimetal and producing a state with electron and hole 
Fermi surfaces, the topological properties of the Weyl semimetal state, i.e. the chiral edge states 
and the corresponding Hall conductivity, survive for moderate I symmetry breaking. Moreover, we 
demonstrate that a new topological phenomenon arises in this case, if an external magnetic field 
along the growth direction of the heterostructure is applied. Namely, this leads to an equilibrium 
dissipationless current, flowing along the direction of the field, whose magnitude is proportional to 
the energy difference between the Weyl nodes and to the magnetic field, with a universal coefficient, 
given by a combination of fundamental constants. 



I. INTRODUCTION 

An interesting recent development in the study of topo- 
logical phases in condensed matter is the realization that 
such phases can exist in gapless semimetalsP^^ This goes 
against the common wisdom that a bulk spectral gap is 
a necessary ingredient in any topological phase, elimi- 
nating the sensitivity to small perturbations that any 
gapless phase would naively exhibit. Such topological 
semimetals are characterized by the presence of two or 
more nondegenerate accidental band-touching points in 
their bandstructure. As first understood by Volovikpffi^ 
these points are characterized by a nontrivial topology, 
and can be thought of as hedgehog-like topological de- 
fects and point sources of Berry flux in momentum space, 
characterized by an integer (±1) topological charge. The 
k ■ p Hamiltonian near such band-touching points takes 
the form of 2 x 2 Hamiltonian of a chiral Weyl fermion, 
hence the name Weyl semimetal. 

In the presence of both time-reversal (TR) and in- 
version (I) symmetries, all bands are doubly degener- 
ate at all momenta in the Brillouin zone (BZ), as fol- 
lows from the Kramers theorem. This means that Weyl 
semimetals require breaking of cither TR or I symme- 
tries,^ which makes a contact between only two bands 
possible, separating the individual Weyl nodes in momen- 
tum space. When Weyl nodes are separated from each 
other in momentum space, they can not be hybridized, 
which makes them indestructible, as they can only dis- 
appear by mutual annihilation of pairs with opposite 
topological charges. This is the mechanism of topologi- 
cal stability of a Weyl semimetal state, which is distinct 
from the spectral-gap protection in insulating topological 
phases 

Perhaps the most interesting variety of Weyl semimet- 



als is obtained when TR is broken. As first shown in 
Ref. HI the simplest kind of a Weyl semimetal can be 
realized in this case, with only two, the minimal num- 
ber required by the fermion-doubling theorem,^ Weyl 
band-touching nodes. This is achieved in a magnetically- 
doped multilayer heterostructure, consisting of thin lay- 
ers of topological insulator (TI) material, with normal 
insulator (NI) spacer layers in between. This particu- 
lar realization of a Weyl semimetal is characterized by a 
semi-quantized anomalous Hall (AH) conductivity, pro- 
portional to the separation between the Weyl nodes in 
momentum space, and chiral surface states, and can be 
thought of as the closest three-dimensional (3D) analog 
of a 2D integer quantum Hall state. We will limit our 
discussion to this system henceforth. 

While the broken TR is necessary to separate the Weyl 
nodes in momentum space and create a topologically- 
stable phase, the remaining I symmetry also plays an 
important role. This symmetry guarantees that the two 
Weyl nodes, separated along the growth direction of the 
TI-NI multilayer in momentum space, occur at the same 
energy. This makes it possible for the Fermi energy to 
coincide with both Weyl nodes simultaneously and thus 
realize a true nodal semimetal phase. The goal of this pa- 
per is to explore the consequences of breaking the inver- 
sion symmetry in the TI-NI multilayer Weyl semimetal, 
in addition to breaking TR. As we demonstrate, a moder- 
ate breaking of inversion symmetry does not destroy the 
topological nature of the state, even though it generally 
moves Weyl nodes to different energies, thus destroying 
the nodal semimetal and producing a state with electron 
and hole Fermi surfaces. However, topological proper- 
ties survive, as long as the I breaking is not too strong. 
Moreover, we demonstrate that a new topological phe- 
nomenon arises in the case, when the Weyl nodes are 
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separated in both energy and momentum. Namely, an 
external magnetic field, applied along growth direction 
of the heterostructure, produces a nondissipative ground 
state current. This current is associated exclusively with 
the zero-mode Landau levels (LLs) of the system and has 
a topological character, in the sense made precise below. 

The rest of paper is organized as follows. In section |TT| 
we present a general Chern-Simons-like description of a 
Weyl semimetal, which demonstrates, that in addition to 
the anomalous Hall effect, there exists another topolog- 
ical effect, namely the magnetic-field driven equilibrium 
current, which arises when the Weyl nodes are separated 
not only in momentum space, but also in energy. In 
section [TTI| we present a model of a TI-NI multilayer het- 
erostructure with broken TR and inversion symmetries 
and discuss the general phase diagram of this system. In 
section [Tv| we present an explicit linear-response calcula- 
tion of the magnetic-field driven equilibrium current in 
the TI-NI multilayer. We conclude in section [V] with a 
discussion of the experimental observability of the pro- 
posed effect. 



Hall conductivity, given by: 
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where 2fco is the separation between the Weyl nodes in 
momentum space. This result can be obtained from the 
Chern-Simons action Eq. ([!]). Indeed the current density 
can be calculated as: 
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where ip = A . Thus in this case 6 — a xy . 

Now let us take 6 M to be purely time-like, i.e. 6 M = 
(6,0,0,0). Taking a functional derivative of the Chern- 
Simons-like action with respect to the vector potential, 
as above, we obtain: 



II. CHERN-SIMONS-LIKE DESCRIPTION OF 
THE WEYL SEMIMETAL 

We will start with a general discussion of the topo- 
logical response in Weyl semimetals, based on 3 + 1- 
dimensional Chern-Simons-like theory, which nicely uni- 
fies bot h the semi-quantized AH response, discussed be- 
forep^ and the new topological magnetic field-induced 
current effect, proposed in this paper. The fact that 
chiral Weyl fermions, separated in momentum and en- 
ergy, induce a Chern-Simons-like term in the action of 
the electromagnetic field, is well-known in the particle- 
physics hteraturep^Hl where it arises in the context of 
CPT and Lorentz-violating extensions of the Standard 
Model. The action, written in the real time representa- 
tion, has the following form: 
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where A^ is the electromagnetic gauge potential and 
(i = 0,1,2,3. Here we use the standard notation for 
the four- vectors x^ — (x°,x), x M = (x ,— x), where x° 
is the temporal component. The four- vector 6 M selects 
a direction in the 4-dimensional space-time, thus violat- 
ing Lorentz invariance. Its temporal component, 6°, is 
proportional to the energy difference between the Weyl 
nodes (we restrict ourselves to the case of only two Weyl 
nodes), while its spatial component, b, is proportional 
to the momentum-space separation between them.EsHini 
Let us first assume that ¥ L has no temporal component. 
By a spatial rotation we can then bring it to the form 
V L = (0,0,0,6), which corresponds to Weyl nodes, sep- 
arated along the z-direction in momentum space. As 
shown in Ref.|H such a state is characterized by a nonzero 
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Taking a — 3, we then obtain; 
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where B = V x A is the magnetic field. The physical 
meaning of this result is that a magnetic field, applied 
to a Weyl semimetal with an energy separation between 
the nodes, will induce a current in the direction of the 
applied field. This current exists in equilibrium and is 
thus nondissipati ve. T his result is known in the particle 
physics literature] 21 * 22 ! in the context of parity-violating 
models of massless chiral particles, as chiral magnetic 
effect. As we will demonstrate below, it is realized in a 
simple model of a Weyl semimetal with broken TR and 
I symmetries. 

A couple of comments are in order here. First, it may 
appear from the above discussion, that energy separation 
between the Weyl nodes is enough in order to obtain the 
dissipationless current and no momentum space separa- 
tion is needed. While this is true in an abstract phe- 
nomenological model, in which Weyl fermions of opposite 
chirality are simply introduced by hand, it is certainly not 
true in our microscopic model, where Weyl fermions ap- 
pear only when separated in momentum space, except at 
a single point in parameter space, corresponding to the 
TI-NI transition in the absence of TR and I breaking. 
Correspondingly, as will be shown below, the time-like 
parameter 6 in the Chern-Simons theory is only nonzero 
when the Weyl nodes are separated in momentum space, 
except if the multilayer system is fine-tuned to the TI- 
NI transition point in the absence of TR-breaking [see 
Eq. (53) below]. Second, Eq. |5| is explicitly isotropic 
in space, i.e. the magnetic field and the corresponding 
dissipationless current can have any direction. Our mi- 
croscopic model, on the other hand, is anisotropic (there 
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is a preferred direction — the growth direction of the mul- 
tilayer) and only becomes effectively isotropic at low en- 
ergies, long distances and weak applied fields. Thus, in 
all our calculations below we assume that the magnetic 
field is applied along the growth direction of the multi- 
layer, since in this case all the calculations can be done 
exactly, without making any additional assumptions of 
weak fields and low energies. The dissipationless cur- 
rent, however, will exist if the magnetic field is applied 
along any other direction, in accordance with Eq. ([5]). 



III. TI-NI MULTILAYER WITH BROKEN 
TIME-REVERSAL AND INVERSION 
SYMMETRIES 

We consider a realization of a Weyl semimetal, based 
on a magnetically-doped TI-NI multilayer heterostruc- 
ture, introduced by one of us in Ref. HJ 
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The first term in Eq.([7]) describes the two (top and bot- 
tom) surface states of an individual TI layer, vp is the 
Fermi velocity, characterizing the surface Dirac fermion, 
which we take to be the same on the top and bottom sur- 
faces of each layer. kj_ is the momentum in the 2D surface 
BZ (we use Ti = 1 units), cr is the triplet of Pauli matri- 
ces, acting on the real spin degree of freedom, and r are 
Pauli matrices, acting on the which surface pseudospin 
degree of freedom. The indices i,j label distinct TI lay- 
ers. The second term describes exchange spin splitting 
of the surface states, which can be induced, for exa mple, 
by doping each TI layer with magnetic impurities! 23 ^ 24 ! 
The remaining terms in Eq.Q describe tunneling be- 
tween top and bottom surfaces within the same TI layer 
(the term, proportional to Ag), and between top and 
bottom surfaces of neighboring TI layers (terms, propor- 
tional to Ad). Longer-range tunneling is assumed to 
be negligible. Partially diagonalizing Eq. Q by Fourier 
transform, we obtain a momentum space Hamiltonian: 



HQs.) = v f t z (z x cr) • k + mo z + A(k z 
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where A(k x ) = A s t x + ±(A D T+e lk * d + H.c). 

The TI-NI multilayer, described by Eq. ([8]), possesses 
inversion symmetry with respect to an inversion center 
placed midway between the top and bottom surfaces in 
any TI or NI layer. Explicitly, the inversion operation 
changes the sign of the momentum and interchanges the 
top and bottom surfaces in each TI layer (spin cr is a 
pseudovector and is thus invariant under inversion): 

l:H(k)^T x H(-k)T x . (9) 

It is easy to see that Eq. ((8|) is indeed invariant under I. 



We now relax the assumption that the multilayer is 
inversion-symmetric, which is perhaps more realistic, as 
some degree of structural inversion asymmetry will be in- 
evitably introduced during growth. Note that inversion- 
asymmetric TI-NI multilayer model (but with preserved 
TR symmetry) was considered before in Ref. |8] 

We thus need to add terms to the multilayer Hamil- 
tonian Eq. ^ , which violate inversion symmetry. There 
are of course infinitely many such terms one can write 
down and we will restrict ourselves only to terms, inde- 
pendent of momentum (i.e. leading-order terms, arising 
in the expansion of a general I-breaking term with respect 
to the momentum). A simple analysis shows that there 
are only two such terms: Vt z and Xt v <j z , which break 
inversion, but respect all other symmetries of the multi- 
layer, in particular rotation around the z-axis. The first 
term corresponds to an electrostatic potential difference 
between the top and bottom surface of a TI layer in each 
unit cell (note that the total electric field in each unit cell 
vanishes). The second term is momentum- independent 
spin-orbit interaction term, which is allowed by the bro- 
ken inversion symmetry. This term is closely analogous 
to the Dresselhaus spin-orbit interaction terms in semi- 
conductors lacking inversion symmetry.^ It is simpler 
than, say, the well-known Dresselhaus terms in semicon- 
ductors with zinc-blende structure due to the fact that 
our system has a uniaxial anisotropy. 

To understand the effect of the I-breaking terms on the 
physical properties of the multilayer, it is useful to start 
from the case when neither TR nor I are violated. In 
this case, diagonalizing Eq. Q gives four bands, which 
are pairwise-degenerate at every momentum, as required 
by the Kramers theorem: 



e ± (k)=±Jv 2 F (k 2 + k 2 y ) + A 2 (k z ), 



(10) 



where A{k z ) = ^A| + A 2 D + 2A S A D cos(fc z d) and d is 
the superlatticc period. Let us now break I without 
breaking TR, i.e. add the Vt z and Xr y a z terms to the 
Hamiltonian, but keep m = 0. In this case the band dis- 
persion can still be found analytically and one obtains: 

4(k) - v 2 F (k 2 x + ky) + V 2 + X 2 + A 2 (k z ) 



± 2\Jv 2 F {V 2 + X 2 )(k 2 x + fc 2 ) + -X 2 A 2 D [1 - cos(2fc 2 d)]. 

(11) 

Thus the double degeneracy at every momentum is lifted, 
as it should be, except at the I and TR-invariant mo- 
menta k x = k v = and k z = 0,ir/d. The role of the 
Vt z term is to split the degeneracy for k x ,k y ^ 0, while 
the role of the Xr v a z term is to split the degeneracy at 
k x = k y = and k z ^ 0, ir/d. In other words, in the ab- 
sence of the Xt v (t z term, the double degeneracy remains 
everywhere along the line k x = k y = in momentum 
space. 

Let us now also break the TR symmetry and add the 
spin-splitting term ma z . In general, it is now impossible 
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to find the energy eigenvalues of the multilayer Hamil- 
tonian analytically, except in a few special cases, which 
we will now discuss, before considering the most general 
situation. The first special case we will consider is A = 0. 
The band dispersion in this case is given by: 

4(k) = v 2 F (k 2 x + k 2 ) + V 2 + m 2 + A 2 (k z ) 

± 2^v 2 F V 2 (k 2 + fc 2 ) + m 2 V 2 + m 2 A 2 (k z ). 

(12) 

The e_ (k) pair of bands exhibits band-touching nodes, 
which are solutions of the equation: 

[v 2 F (kl + K) + m " + A2 (M ~V 2 ] 2 = A(m 2 -V 2 )A 2 (k z ). 

(13) 

Clearly, two cases must be considered: m > V and m < 
V. In the latter case, a solution exists only when A(fc z ) = 
0, which happens when As/ Ad = ±1 at k z = w/d,Q 
correspondingly. Eq. (fl~3|) simplifies to: 



It is obvious already from the form of Eq. ( 19 1 that the 
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This describes a circular nodal line in the icy-plane, 
which, however, exists only at the TI-NI critical point 
Ag /Ad = ±1 (a stable Weyl semimetal phase can be 
obtained in this case, if one includes momentum depen- 
dence of the Agjj tunneling amplitudes, breaking the 
z-axis rotational symmetry; we will not pursue this direc- 
tion here, as it has been investigated in detail in Ref. |8|) . 
In the m > V case, Eq. (13) simplifies to: 



v 2 F (kl + kl) + [v 7 ™ 2 ~ V 2 - A(k z ) 



= 0. 



The solution is: 



k x = k y = 0, A{k z ) = yjm 2 - V 2 . 



(15) 



(16) 



This corresponds to the point-node Weyl semimetal, con- 
sidered in Ref. 01 The effect of the Vt z term in this case 
is only to reduce the distance between the Weyl nodes 
in momentum space. Henceforth, we will assume that 
m > V and, keeping the above analysis in mind, neglect 
the Vt z term entirely. 

Now let us analyze the effect of the I-breaking spin- 
orbit interaction term Xr y a z in the presence of the broken 
TR symmetry due to the raa z term. In this case, the 
band dispersion can not be found analytically in general, 



except at the k x 
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point. At this point we obtain: 



e s± (k z ) = ms ± v/A 2 + A 2 (k z ) - 2s\A D sin(k z d), (17) 

where s = ±. The Weyl nodes are given by the solution 
of the equation: 

e + _(fc z ) = e_+(fc z ), (18) 

or, explicitly: 

2m = v/A 2 + A 2 (fc z ) - 2\A D sin(k z d) 
+ v/A 2 + A 2 (fc z ) + 2AA D sin(fc z rf). (19) 



nodes no longer occur at zero energy. Eq. ( 19 ) can not be 



solved analytically in its general form, but we can solve it 
perturbatively in A. To first order in A, the position of the 
nodes along the z-axis in momentum space in unchanged 
and is given by the equation: 



m = A(fc z 



(20) 



Assuming for concreteness that the tunneling amplitudes 
As : d are both positive, the nodes are then located at 
kf —ir/d± k 0l whereP 
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The shift of the Weyl nodes from the zero energy can 
then be obtained from Eq. (17) by expanding it to first 
order in A. We obtain: 

XAjj sin(k z d) 



Ae ± = 



^{m 2 c2 -m 2 ){m 2 -m 2 cl ), (22) 
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where m c \ — \ A$ — Ad| and m C 2 = As + Ajj are, corre- 
spondingly, the lower and upper critical values of the spin 
splitting r7i, between which the Weyl semimetal phase ex- 
ists 3 

Thus, in the presence of the \r y a z term, the two Weyl 
nodes are shifted in the opposite directions in energy, 
as seen in Fig. [T] It is no longer possible to align the 
Fermi energy with both nodes simultaneously and the 
semimetal at charge neutrality will have equal-volume 
electron and hole Fermi surfaces. The Weyl nodes ex- 
ist as long as the electron and hole Fermi surfaces are 
separated. Upon increasing the magnitude of A, the 
momentum-space separation between the nodes and be- 
tween the electron and hole Fermi surf aces decre ases. 
Eventually, at a critical value A c i = \Jm 2 — m 2 1: the 
Fermi surfaces touch and the Weyl nodes disappear, see 
Fig. [2j The Fermi surfaces still exist, however, until, at 
a still larger value of A C 2 = \J '(m + Ad) 2 — A|, a fully 
gapped state is produced, see Fig. [3] In fact, the phase 
diagram of the system in the m — A plane is somewhat 
more involved, and we will explore it in more detail be- 
low. 

An important question, which arises in this con- 
text, is whether the topological properties of the Weyl 
semimetal, like the chiral edge states and the anomalous 
Hall conductivity, proportional to the momentum-space 
separation between the Weyl nodes, survive when the 
nodes move away from the zero energy and the electron 
and hole Fermi surfaces appear. We will demonstrate be- 
low that this is indeed the case: the system retains its 
topological properties as long as it is gapless, in fact even 
after the Weyl nodes disappear above A = A c i . 

To see this we start from the momentum-space Hamil- 
tonian of the multilayer, which, after a canonical trans- 
formation: 
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FIG. 1. (Color online) A plot of the band dispersion along the 
k x = k y = line for a moderate magnitude of the inversion 
symmetry breaking, corresponding to A/As = 0.3. The Weyl 
band-touching nodes are visibly shifted in opposite directions 
in energy, but are otherwise intact. At charge neutrality the 
Fermi energy is pinned at e F = 0, and there are equal- volume 
electron and hole Fermi surfaces. 




FIG. 2. (Color online) Band dispersion along the k 
line for A 



ky = 



Aci = \/m 2 — (As — Ad) 2 . The electron and 
hole Fermi surfaces touch and the Weyl nodes disappear. 



takes the following form: 

ft(k) = v F {z x a) • k + [m + A{k z )]a z + Xr y . (24) 

We now add a small external magnetic field B along the 
z-direction (i.e. the growth direction of the multilayer), 
which we will set back to zero at the end. The Hamilto- 
nian becomes: 

U(k z ) = v F (zxa)-(-iV + -A)+[m, + K{k z )]a z + \T y , 



where A = xBy in the Landau gauge. To find the Lan- 
dau levels, we introduce LL ladder operators in terms of 
the components of the kinetic momentum it — — iV + -A 
in the standard way (assuming B > for concreteness) : 



(tt x -ITTy), a< 



(n x +m y ), 



(26) 
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where Is — y/c/eB is the magnetic length. In terms of 



FIG. 3. (Color online) Band dispersion along the k, 
line for A = A C 2 = ^(m + Ac 



k LI 







A 2 



The Fermi surfaces 
disappear entirely and the semimetallic state gives way to a 
fully gapped insulator. 



the ladder operators, the Hamiltonian takes the form: 

H(k z ) = ^ (a+a-a-a^) + [m + A(k z )}a z + XTy, (27) 
v2 



where lob — Vf/^b- It is clear from Eq. (27), that its 
eigenstates have the following general form: 

\n) = z n+t \n- l,+,t) + Zn+±\n,+,l) +z n -f\n- l,-,f) 
+ z n ^\n,-,i), (28) 

where |n, ±,tl) is the nth LL eigenstate on the top (+) 
or bottom (-) surface with spin up or down, and z are 
complex amplitudes. We have omitted all the additional 
eigenstate labels on the left hand side for brevity, leaving 
only the LL index n. In the presence of the Xt v term, the 
above Hamiltonian can not be diagonalized analytically 
in general. However, we can find the n = LL dispersions 
analytically, and this is in fact all that is needed for our 
purposes. Indeed, from Eq. (28) it is immediately clear 
that the n — pair of LLs are polarized downwards, i.e. 
we can replace a z = — 1 and then the n = LL dispersion 
is given by the eigenvalues of: 



n n=0 {k z ) = -m - A(k z ) + \t» . 



Diagonalizing Eq. (29), we then find the n 
persions: 



(29) 
LL dis- 



-in 



± VX 2 + A 2 (k z ) + 2\A D sin{k z d). (30) 



The eo-(fcz) level is below the Fermi energy ep — for all 
values of k z . The eo+(k z ) LL, on the other hand, crosses 
the Fermi energy twice at the following momenta: 



i.± * + 1 

kf = — ± -arccos 

d d 



2A D {A% + X 2 ) 



(31) 



where 

f ± = A S (A 2 

±x x /[aJ 



A 2 D + X 2 - m 2 ) 



X 2 - (to- A D ) 2 ][(m + A D ) 2 



A 2 



-A 2 ]. 
(32) 
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This reduces to the already known result kf = n/dztko^ 
where ko is given by Eq. (21 1 in the limit A = 0. As A 
increases, k~ becomes equal to ir/d at the lower-critical 
value A c i = -\/m 2 — (As — A^) 2 , at which the electron 
and hole Fermi surfaces touch. At the upper critical value 
A c2 = \J (m + A D ) 2 — A|, the Fermi surfaces disappear 
entirely, giving rise to a fully gapped state, and the n = 
LL no longer crosses the Fermi energy. 

There is a subtlety, associated with the A C 2(m) critical 
line. Two cases need to be considered separately. When 
As > Ad, i.e. when the multilayer in the presence of 
both TR and I is an ordinary insulator, the X C 2(m) line 
starts at m = m cl and separates a semimetal from an 
ordinary insulator, as shown in Fig. [4] When A^> > As, 
however, the system is a topological insulator when TR 
and I are unbroken. If TR is preserved, but I is broken 
by increasing the magnitude of A, there is e ventually a 
TI to NI transition at a critical value of A = J A 2 -, — A|. 
Correspondingly, the \ c i{m) line has a different shape, 
shown in Fig. [5] 

Finally, a third critical line A C 3 = \J{m — A_o) 2 — A| 
exists, separating 3D quantum anomalous Hall (QAH) 
insulator^ from the semimetallic phase with Weyl nodes 
and Fermi surfaces. 

The most important consequence of the above analysis 
is that the anomalous Hall conductivity: 
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FIG. 4. Phase diagram of the multilayer structure in the m— A 
plane for As > Ac. There are five distinct phases: (1) In- 
sulator (INS), outside of the A C 2 line. (2) Nodal semimetal 
(NSM), which exists along the interval A = 0, m c i < m < 
m C 2. (3) Semimetal with electron and hole pockets, but 
with Weyl nodes preserved (FSSM 1). This is bounded by 
the A c i = ± \/m 2 



and A c3 = ±^/(m — A D ) 2 



A 2 

lines. (4) Semimetal without Weyl nodes (FSSM 2). This 
is bounded by the A C 2 = [m + A_d) 2 — A| and the A c i 
lines. (5) QAH insulator phase, enclosed by the A C 3 line. The 
anomalous Hall conductivity is nonzero in all semimetallic 
phases and the QAH phase. 



is nonzero beyond A = A c i, when the Weyl nodes dis- 
appear, persisting all the way to the A = A C 2 point, 
when a fully gapped insulating state emerges (chiral sur- 
face states disappear when the Weyl nodes disappear at 
A = Ai). This follows simply from the fact that one of 
the two n = LLs always dips below the Fermi energy 
in the interval k~ < k z < kf (the second n = LL 
is always below the Fermi energy) in all gapless phases, 
thus giving rise to a nonzero Hall conductivity (relative 
to the insulating state at m < rn c i, where there is one 
zero-mode LL above and one zero-mode LL below the 
Fermi energy and the Hall conductivity is zero)P The 
connection between the nonzero anomalous Hall conduc- 
tivity and the zero-mode LL crossing of the Fermi energy 
can be seen most easily by referring to the Stfeda formula 
for the quan tized nonclassical part of the anomalous Hall 
conductivity pM3 



DSN 



dB 



(34) 



B=0 



where SN is the extra charge density, corresponding to 
the filled zero- mode LL, when k~ < k z < kf: 



SN 



2^1 j k - 



kt dk r . 



2tt 



e\kt-K)B 
2ixhc 



(35) 



The nonquantized part of the anomalous Hall conduc- 
tivity, on the other hand, is always zero due to charge 



neutrality as long as the Fermi level remains pinned at 
ep = 0. The full phase diagram of the multilayer in the 
m — A plane is shown in Figs. HI [5] 



IV. MAGNETIC FIELD DRIVEN 
EQUILIBRIUM CURRENT 

In this section we will demonstrate that, in addition to 
the semi-quantized AH conductivity and chiral "Fermi- 
arc" surface state, the TI-NI multilayer with broken TR 
and inversion is characterized by another topological phe- 
nomenon, namely the external magnetic field driven per- 
sistent current, which we have introduced in section [H] 
using general arguments, based on Chern-Simons like de- 
scription of the Weyl semimetal. Here we will limit our 
analysis to small values of A, in which case all calcula- 
tions can be done analytically, and which is perhaps the 
most relevant regime experimentally. 

We will calculate the equilibrium current in linear re- 
sponse with respect to the Xr v a z term. In this case, the 
energy separation between the Weyl modes: 



Ae = |Ae_ 



Ae+| 



A 



A s |m| 



\] ( TO c2 — m 2 )(m 2 — m 2 x ), 



(36) 

can be thought of as an "electrochemical potential dif- 
ference" between them. We thus want to calculate an 
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QAH 



FIG. 5. Same as in Fig. [4] but for A D > As- There is a ad- 
ditional transition in this case, which occurs as A is increased 
along the m = line. The |A| < A C 2 1 phase is a TI, while the 
|A| > | A e 2 1 phase is an ordinary insulator. 



electrical current that arises in response to this "electro- 
chemical potential difference" . Note that we have denned 
Ae as the energy shift of the node of positive chirality 
(which is at k~ for m > but at fc+ for m < 0) minus 
the energy shift of the node of negative chirality: defined 
this way, Ae is independent of the sign of m. We also 
assume that an external magnetic field of magnitude B 
is applied along the growth direction of the multilayer, 
which we take to be the z-direction (note that the spin 
magnetization m does not contribute to the total internal 
orbital magnetic field for a sample of an infinite slab ge- 
ometry, due to perfect cancellation with the correspond- 
ing demagnetizing field) . We will ignore the Zccman con- 
tribution of the external field, assuming it to give only 
a small correction to m. The unperturbed Hamiltonian, 
after the canonical transformation of Eq. (23), is given 
then by: 



n(k z ) =v F (zxcr)-l-iV 



-A(k z )]a*. (37) 



Diagonalizing A(fc z ), which is a constant of motion, 
brings the Hamiltonian to a block-diagonal form, with 
two independent 2x2 blocks: 



H±(k z ) = v F (z x <x) ■ (-iV 



m±(k z )a\ (38) 



where m±(k z ) — m ± A(fc 2 ). Using Landau gauge A = 
xBy (for convenience, we will implicitly assume B > 
in all calculations, but the final result does depend on 
the sign of B) and introducing LL ladder operators as in 
section IIIIl we obtain: 



T~L a (k z 



la' a — a a' 



m a {k z )a z 



(39) 



where a = ±. It is clear that the eigenstates of each 2x2 
block of Eq. ( 39 ) have the following general form: 



lA)+v$ sa \n,l). 



(40) 



Here n = 0,1,2,... is the LL index and ,s = ± labels 
the electron-like and hole-like LLs correspondingly. The 
eigenvalues of Eq. ( 39 ) are easily found and are given by: 



£ nsa (k z ) = sJ2uj%n + m%(k z ), n > 1, 

e 0a = ~m a (k z ). (41) 
The corresponding eigenvectors are: 



n > 1, 
(42) 

We note the following symmetry properties of the eigen- 
functions Eq. ( 42 ) , which will play an important role in 
our analysis: 




J n I 



(v+ a '\v: 



while: 



n > 1, 



n > 1. 



(43) 



(44) 



The total four-component spinor cigenfunctions of 
Eq. (39), which we will denote simply by \nsa), are given 
by the tensor product of \v^ a ) and the eigenf unctions of 
the A(k z ) operator: 



i.e. 



I a \ 1 (, A s + A D e-^ d 



|n*or) = \u a ) <g> \v%*). 



(45) 



(46) 



Let us point out the following matrix elements of the 
eigenstates of the tunneling operator, which will be im- 
portant in the calculation below: 



{u + \t x \u-)=z 



Ad sia(k z d) 
A(fc z ) ' 



i +i v\ -\ .A s + A D cos(fc z d) 

)=t ^ . (47) 

We evaluate the electrical current in the z-direction in 
linear response with respect to the inversion-symmetry 
violating term Xr v a z . The standard expression is: 

A 



Jz 



2*e B L z 



E 



(nsa\j z \n s a 



n,n's,s'a,a'k z 



x (n s a \T v \nsa) -, (48) 

where L z is the size of the system in the growth direction, 
we have suppressed explicit dependence on k z everywhere 



for brevity, and we are using the momentum space Hamil- 
tonian Eq. (37), i.e. after the canonical transformation 



of Eq. ( 23 1 . The current operator j z is given by: 



&H{k z 

' dk z 



eA D da z [t x sm(k z d) + t v cos(k z d)} . 

(49) 

To evaluate Eq. (48), we first note that both j z and t v 



operators are diagonal in the spin indices, which means, 
taking into account Eq. ( 40 ) , that the matrix elements of 



these operators, appearing in Eq. ( 48 ) , are diagonal in the 
LL index n. Furthermore, if the Fermi level is pinned by 
charge neutrality at zero energy ep = 0, only the terms 



with s ^ s' contribute to Eq. (48) due to the difference 



of the Fermi distribution functions in the numerator (we 
assume the temperature to be lower than all other energy 
scales in the problem). Then we obtain: 



A 



Jz 



E 



2nl%L z — 

n,S7=s ,a,at ,K Z 



[nsa\j z \ns' a' 



x (ns a \T y \nsa) - - -, (50) 



To simplify this further we use Eqs. (43 1, (44) and (47) 



According to these, the product of the matrix elements of 



j z and t v , appearing in Eq. (50), is antisymmetric with 



respect to the interchange of the a and a' indices for all 
n > 1. This means that the contribution of all n > 1 LLs 
cancels out and only the two n — LLs contribute to the 
expectation value of the current. Then we obtain: 

3*=2=FT E (0a\j z \0a'){0a'\ry\0a) 

B z a^a'.k z 

^F(eoa) - n F (e 0a ') 

X 

^0a ~~ ^Oa' 



(51) 



Since eo+(fc z ) = —m + {k z ) = —m — A(k z ) is always neg- 
ative, while eo_(fc z ) = — to_(/c 2 ) = —to + A(fc 2 ) changes 
sign at the locations of the Weyl nodes, only the subset 
of the ID BZ, in which m + (k z ) and TO_(fc z ) have oppo- 
site signs, contributes to Eq. (51 1 . Evaluating the matrix 
elements in Eq. (51 1, we then obtain: 



Jz 



eA D d r/ d+k ° dk z 
2n£% J v/d _ ko ~2^ 



A 2 {k z ) cos(k z d) + A S A D sin 2 (k z d) 
X A3(fc z ) ' [b2) 

where we have used the fact that the above integral, eval- 
uated over the whole BZ, vanishes. Evaluating the inte- 
gral over the interval between the locations of the Weyl 
nodes in Eq. (53), we obtain: 



Jz 



4ir 2 \m\Asc 



\f( m c2 ~ m 2 )(m 2 — to^). (53) 



Using Eq. ( 36 ) for the energy difference between the Weyl 
nodes, induced by the inversion-breaking term Xr y a z , we 



finally obtain: 



Jz 



e 2 Ae 
4ir 2 h 2 c 



B, 



(54) 



where we have restored explicit fi. Note that B in 
Eq. (54) can have both signs, and the current correctly 
changes sign under TR, since Ae is TR-invariant. There 
is thus indeed a purely equilibrium nondissipative cur- 
rent, as predicted by the Chern-Simons theory analysis of 
section||l] The current is proportional to both the applied 
external field and to energy separation (for small separa- 
tion) between the Weyl nodes. The remaining coefficient 
is a universal combination of fundamental constants. 

The above result can be understood physically 
by appealing to the well-known Adler-Bell-Jackiw 
anomalyP^ni Following the argument by Nielsen and Ni- 
nomiyap2122l j e t us imagine that there is an electric field 
E, applied in the same direction as the magnetic field B. 
In the presence of the magnetic field, the electric field will 
lead to transfer of particles between the two Weyl nodes 
through the zero-mode LL, which has a definite chirality 
near each of the Weyl nodes, reflecting the chirality of 
the nodes. The rate of the particle transfer between the 
left ( — ) and right (+) node per unit volume is given byP^ 



d(A_ - N+) eE 1 



dt 



2kTi 2ir£ 2 2-n 2 h 2 . 



-EB. 



(55) 



Since the Weyl nodes have an energy difference Ae, this 
particle transfer process has an associated power cost per 
unit volume: 



P = A( 



d{N_ -N+)/2 



(56) 



This power is provided by the current of magnitude j, so 
that P = jE. Thus we obtain: 



jE = Ae 



ATT 2 h 2 C 



EB. 



(57) 



Cancelling E on both sides of this equation and sending 
it to zero, we obtain: 



7 = o~B, 

4n 2 h 2 c 



(58) 



which is identical to Eq. ( 54 ) (the sign of the current is 



not determined by the above energy balance argument) 



V. DISCUSSION AND CONCLUSIONS 

Let us now address the experimental observability 
of the proposed magnetic-field-driven persistent current. 
First let us estimate possible current density magnitudes. 
We take the energy separation between the Weyl nodes 
Ae ~ 1 meV and the magnetic field B~ 1 T. This gives 
j z ~ 0.1 A/cm . This is certainly an easily measurable 
current. 
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The effect itself will of course exist only at low enough 
temperatures and in clean enough samples, so that, in 
particular, the temperature and the impurity scattering 
rate are smaller than Ae, the energy separation between 
the Weyl nodes. This implies temperatures less than 
about 10 K, if we take Ae ~ 1 mcV. Thus the proposed 
effect may be observable. The biggest unknown here is of 
course the value of the spin-orbit interaction parameter 
A. In the above estimate we have assumed it to be of 
order ~ 1 meV. It may, however, turn out to be much 
smaller, in which case the effect we have described may 
turn out to be unobservable. 

In conclusion, we have studied Weyl semimctal, real- 
ized in a TI-NI multilayer heterostructure, in which both 
TR and I symmetries are violated. Breaking TR symme- 
try moves the Weyl nodes to different points in momen- 
tum space, while breaking inversion shifts them to differ- 
ent energies, thus producing a state with equal-volume 
electron and hole pockets at charge neutrality. We have 



demonstrated that if an external magnetic held is applied 
to this system along the direction of the vector, connect- 
ing the Weyl nodes, an equilibrium nondissipative cur- 
rent flows in this system. The current is proportional to 
both the applied magnetic field and the energy difference 
between the Weyl nodes, with a coefficient that is a uni- 
versal combination of fundamental constants. This cur- 
rent may be measurable under experimentally-achievable 
conditions. In a sample, not connected to external leads, 
a voltage in response to an applied magnetic field will in- 
stead be measured. This effect can potentially be useful 
for magnetic field detection. 
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